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Symmetry-breaking magnetization dynamics of a spin-1 Bose-Einstein condensate (BEC) due to 
the dipole-dipole interaction are investigated using the mean-field and Bogoliubov theories. When a 
magnetic field is applied along the symmetry axis of a pancake-shaped BEC in the m — hyperfine 
sublevel, transverse magnetization develops breaking the chiral or axial symmetry. A variety of 
magnetization patterns are formed depending on the strength of the applied magnetic field. The 
proposed phenomena can be observed in 87 Rb and 23 Na condensates. 
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I. INTRODUCTION 

A Bose-Einstein condensate (BEC) of atoms with spin 
degrees of freedom (spinor BEC) allows the study of mag- 
netism in a quantum fluid. There are two mechanisms 
of magnetization for a spinor BEC: the ferromagnetic 
contact interaction and the magnetic dipole interaction 
(MDI) between atoms. Magnetization dynamics due to 
the ferromagnetic contact interaction have been observed 
for a spin-1 87 Rb BEC However, magnetization dy- 

namics due to the MDI have not been studied yet, and 
this is the subject of the present paper. 

A BEC of 52 Cr atoms with a large magnetic dipole mo- 
ment (6 fiB with /iB being the Bohr magneton) has been 
realized by the Stuttgart group [3[ and its anisotropic 
behaviors originating from the anisotropy of the dipole- 
dipole interaction have been observed 0, @, 0, Q. 
While the magnetic dipole moment of spin-1 alkali atoms 
(/zb/2) is much smaller than that of 52 Cr, it is neverthe- 
less predicted that a small MDI can create spin textures 
in a 87 Rb BEC [1, [IqI El • The crystalline magnetic order 
observed in a 87 Rb BEC is considered to be caused by 
the MDI 111- MDI effects have been detected in 39 K [l3[ 
and 7 Li [14[ BECs using Feshbach resonance. 

In this paper, we show that magnetization dynami- 
cally develops in spin-1 87 Rb and 23 Na BECs due to the 
MDI. We consider a situation in which a BEC in the 
m = magnetic sublevel is confined in an axisymmet- 
ric pancake-shaped trap and a magnetic field is applied 
along the symmetry axis. We numerically solve the time- 
dependent nonlocal Gross-Pitaevskii (GP) equation in- 
cluding the MDI and show that magnetization develops 
in the direction perpendicular to the magnetic field due 
to the MDI. These magnetization dynamics break the 
chiral or axial symmetry spontaneously. We find that 
various magnetization patterns emerge depending on the 
strength of the magnetic field. For spin-1 87 Rb, we can 
suppress magnetization due to the ferromagnetic contact 
interaction using the microwave-induced quadratic Zee- 
man effect [15|, and the pure MDI effect can thus be 
observed. We perform a Bogoliubov analysis and show 
that the magnetization is triggered by the dynamical in- 



stability. We also employ the variational method with 
the Gaussian approximation to explain the numerical re- 
sults. 

This paper is organized as follows. Section U formu- 
lates the mean-field and Bogoliubov theories to study 
the present system. Section Hill numerically studies the 
Bogoliubov spectra and demonstrates the magnetization 
dynamics. Section HVl analyzes the phenomena using the 
variational method. Section [VI gives conclusions to the 
study. 



II. FORMULATION OF THE PROBLEM 

A. Mean-field theory 

We consider a system of spin-1 bosonic atoms with 
mass M and magnetic dipole moment confined in 

an axisymmetric harmonic potential V(r) = M{uJ\r\ + 
u 2 z 2 )/2 with r± = (x 2 + y 2 ) 1 ^ 2 . A uniform magnetic 
field B z is applied in the z direction and the linear Zee- 
man energy is given by —jiBB z m/2 for magnetic sub- 
levels m = 0, ±1. We neglect the magnetic quadratic 
Zeeman effect, since the strength of the magnetic field 
considered here is B z < 10 mG. Instead, we assume that 
the microwave- induced quadratic Zeeman effect [15] lifts 
the energy of the m = ±1 states by Q. 

We employ the mean-field approximation and the con- 
densate is described by the macroscopic wave function 
-0 m with magnetic sublevel m, which satisfies the nor- 
malization condition J |?/; m | 2 <ir = TV, with N being 
the total number of atoms. The nonlocal GP equations 
including the MDI are given by 
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Where p = Em l^m| 2 , F = £ mm , 4>m(f)rnm>tl>m> , / is 

the vector of the spin-1 matrices, F± = i 7 ^ ± ii 7 ^, and 
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The spin-independent and spin-dependent contact- 
interaction parameters go and g\ have the forms 
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where as is the s- wave scattering length for colliding 
atoms with total spin S. The MDI produces an effec- 
tive magnetic field: 
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where //q is the magnetic permeability of vacuum and 
e = (r — r f )/\r — r f \. Equation ([1]) is numerically solved 
in Fourier space for the kinetic term and in real space 
for the other terms using a fast Fourier transform (FFT). 
The FFT is also used to calculate the convolution integral 
in Eq. g]). 

In the present paper, we consider an initial state in 
which all the atoms are in the m = sublevel. In order 
to simulate this situation, we prepare the ground state 
of ipo with ip±i = by the imaginary-time evolution of 
Eq. (pQ). Small noise (a complex random number on each 
mesh) is then applied to the initial state of ip±i to break 
the symmetry and trigger magnetization due to dynam- 
ical instability. The small noise corresponds to quantum 
fluctuation, thermal atoms, and residual atoms in an ex- 
periment. We note that if the initial state of m = ±1 



is i/j±i = 0, the right-hand side of Eq. (|lb|) vanishes and 
ip±i never develops within the mean-field theory. Mag- 
netization thus occurs if small noise in the m = ±1 state 
is exponentially amplified by dynamical instabilities. 



B. Bogoliubov analysis 

We investigate the Bogoliubov excitation spectrum of 
magnons for the stationary state of ipo. Assuming that 
ip±i is small and neglecting the second and higher orders 
of ip±i in Eq. (pQ), we obtain 
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where = /io/i|/(167r) and e± = e x ± ie y . Using the 
mode functions u±\ and v±i, we write a single- mode ex- 
citation of i/j±i as 
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where <p = arg(x + iy), L is an integer, and \i is the 
chemical potential, 
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Each ip±i in Eq. ([6]) has 2(L — l)-fold symmetry 
around the z axis. Substituting Eq. (|6j) and iJjq = 
\i/jq\ exp(— ifit/h) into Eq. ([5]) yields the closed form of 
the nonlocal Bogoliubov-de Gennes equations, 



2M 



Vi 



2 (^±1-1) 21 



V T P + Q + <?o|^o| 2 - fJL U±i + 91 |^o| 2 (u±i + v Tl ) 



\r — r' 



2M 



Vi 



2 (L T l-lf 



V T P + Q + 9o\H -li>v ±1 + gi\M («=fi + v ±i) 



"f / dr 'y^f ~ 3e ') Mr')+v ±1 (r')}+3el [ u±1 (r') + v Tl (r')}} |^(r)| = -tiuv ±1 , 



(8a) 



(8b) 



where V\ z = d% ± +r A _ 1 d r± + d 2 z . We numerically diag- onal bases, e.g., the harmonic-oscillator eigenfunctions. 
onalize Eq. J5]) by expanding u±\ and v±\ with orthog- 
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FIG. 1: (color) Imaginary part of Bogoliubov frequency oo 
for magnon excitation from the m = state of 87 Rb atoms as 
a function of P — Q with Q = 50hu± . The excitation mode 
has 2(L — l)-fold symmetry around the z axis, where L is 
defined in Eq. ©. The number of atoms is N = 10 5 and the 
radial and axial trap frequencies are co± = 2tv x 100 Hz and 
u z = 2tt x 400 Hz. The values of P indicated by the arrows 
are used in Figs. [2] and [3] 

If complex frequencies uo emerge, the stationary state 
becomes dynamically unstable against excitations of 
magnons. 

III. MAGNETIZATION DYNAMICS AND 
BOGOLIUBOV SPECTRA FOR ALKALI ATOMS 

A. Spin-1 rubidium 87 

We first consider a spin- 1 87 Rb BEC. The scattering; 
lengths are ao = 101. 8<2b and a2 = 100. 4<2b [16] with 
<2b being the Bohr radius. The spin-dependent contact- 
interaction parameter g\ is therefore negative and the 
ground state is ferromagnetic [1]. In order to distinguish 
magnetization by the MDI from that by the ferromag- 
netic contact interaction, we apply microwave radiation 
to lift the energy of the m = ±1 states by Q, which must 
be much larger than the ferromagnetic interaction energy 
\gi\p. The magnetization due to the ferromagnetic con- 
tact interaction is thus suppressed and the pure effect of 
the MDI can be observed. 

To investigate the dynamical instability against mag- 
netization, we numerically diagonalize Eq. (|SJ). Figure [1] 
shows the imaginary part of the Bogoliubov frequencies 
uo as a function of the applied magnetic field, where 
N = 10 5 atoms are confined in a pancake-shaped trap 
with uj ± = 2tt x 100 Hz and u z = 2tt x 400 Hz. The 
microwave- induced quadratic Zeeman energy Q is cho- 
sen to be 50hujj_, which is sufficient to suppress mag- 
netization by the ferromagnetic contact interaction. In 
fact, we have confirmed that the Bogoliubov spectrum 
is always real if the MDI is absent, = 0, for this 
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FIG. 2: (color) Time evolution of the population of the m — 
1 state / \i/j!\ 2 dr/N for the values of p = (P - Q)/(huj±) 
indicated by the arrows in Fig. [T] The parameters are the 
same as those in Fig.[T] The populations of the m = — 1 state 
are too small to be discerned at this scale of the ordinate. 



value of Q. The linear Zeeman energy P = 50huj± cor- 
responds to B z ~ 7.15 mG. From Fig. [TJ we find that 
magnon modes with various rotational symmetries (var- 
ious L) become dynamically unstable, depending on the 
linear Zeeman energy P. There is no imaginary part for 
p = (P — Q)/(huj±) < —0.1, while many peaks in the 
imaginary part exist for p > 0.5. 

Figure [2] shows time evolution of f \^i\ 2 dr/N for the 
values of P indicated by the arrows in Fig. [TJ The tran- 
sition from the m = state to the m = 1 state occurs 
due to the dynamical instability shown in Fig. [TJ From 
Fig. [2J we find that the transition occurs periodically ex- 
cept for p = 0.42 (green line). The complicated behavior 
for p = 0.42 originates from the fact that the dynami- 
cally unstable modes are not only L = 2 but also L = 3 
(see Fig. [1]). We note that the transition to the m = — 1 
state is negligibly small and the total spin in the z direc- 
tion J(\i/ji\ 2 — |^_i| 2 )<ir is not conserved, indicating that 
the transition is not due to the spin-exchange contact 
interaction but due to the MDI. Since the z component 
of the total angular momentum must be conserved, the 
system acquires orbital angular momentum. The trans- 
fer of the spin angular momentum to the orbital angular 
momentum in a spinor dipolar BEC also occurs in the 
Einstein-de Haas effect [13, El, EH - 

Figure [3] shows transverse magnetization at the times 
of the first peaks of J \ipi\ 2 dr (the first peaks of the lines 
in Fig. [2]) for the linear Zeeman energies P indicated by 
the arrows in Fig. [TJ A variety of magnetization patterns 
with 2(L — l)-fold symmetry emerge depending on the 
strength of the applied magnetic field. The closure struc- 
ture of the magnetization in Fig. [2] (a) is an energetically 
favorable structure for the MDI energy. The directions of 
the magnetization vectors in the closure structure have 
clockwise and counterclockwise symmetry, and therefore 
the spin- vortex generation in Fig. [3] (a) breaks the chiral 
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FIG. 4: (color) Imaginary part of the Bogoliubov frequency 
uj for magnon excitation from the m — state of 23 Na atoms 
as a function of the linear Zeeman energy P. The range of 
P is (a) < P/{hujA_) < 0.5 and (b) 1.5 < P/(Huj ± ) < 2. 
The excitation mode has 2(L — l)-fold symmetry around the 
z axis, where L is defined in Eq. JB}. The number of atoms 
is N = 10 6 and the trap frequencies are the same as those in 
Fig. [1] The values of P indicated by the arrows in (b) are 
used in Fig. [5] 



symmetry. The m = 1 component of these spin vor- 
tices is ipi oc e~ %< ^ . This situation is different from the 
spin-vortex generation by the ferromagnetic contact in- 
teraction, in which polar-core vortices of ip±i oc e ±l ^ and 
ip±i ex emerge with an equal probability [20]. The 
closure structures are also seen in Figs. [3] (b)-|3] (e). The 
magnetization in Figs. [3] (b)-|3] (e) caused by the dynam- 
ical instability with L > 1 exhibits a variety of patterns, 
breaking the axisymmetry of the system. 
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FIG. 3: (color) Magnitude of the integrated transverse 
magnetization | J F+dz\ (right panels) and its direction arg 
(J F+dz) (left panels) at the first peaks of the curves in Fig. [2] 
The values of p = (P — Q)/(huj±) used are indicated by the 
arrows in Fig.E The unit of / F + dz is NMu ± /h. The length 
of the vector is proportional to | f F+dz\. The field of view is 
9.7 x 9.7 fim. 



B. Spin-1 sodium 23 

Next we consider a spin-1 23 Na BEC. The scattering 
lengths are given by (ao + 2*22) /3 = 53.4<2b [21] and 
&2 — &o = 2.47aB [22j]. The spin-dependent contact- 
interaction parameter g\ is then positive and the polar 
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(a) P/h(0j=1 .63, <B L t=640 (b) p/hco ± =1 .875, co>670 
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FIG. 5: (color) Integrated transverse magnetization | J F+dz\ 
at the time when f\ipi\ 2 dr becomes the first maximum in 
time evolution for the values of P indicated by the arrows in 
Fig. H (b). The unit of | f F+dz\ is NMcv ± /h. The field of 
view is 21 x 21 /mi. 



state (m = 0) is energetically favorable. Spontaneous 
magnetization due to the contact interaction is therefore 
suppressed and the microwave-induced Zeeman effect is 
unnecessary (Q = 0). The number of atoms is assumed 
to be N = 10 6 and the trap frequencies are the same as 
those in Sec. IIII Al 

Figure 2] shows the imaginary part of the Bogoliubov 
frequency obtained by numerically diagonalizing Eq. ([8]). 
Compared with the case of 87 Rb in Fig. [TJ the width and 
height of the peaks are small for < P/(Tiuj±) < 0.5 
[Fig. [31 (a)]. The width and height of the main peaks 
gradually increase and saturate for P/(hu±) ~ 2 [Fig. [H 
(b)]. For P < 0, there is no imaginary part. 

We numerically solve the GP equation for the values of 
P indicated by the arrows in Fig. [H (b). The initial state 
is prepared by the same method as for 87 Rb. Figure [5] 
shows the integrated transverse magnetization | J F+dz\ 
at the time of the first peak of J \ipi\ 2 dr in the time 
evolution. Many radial nodes in the magnetization pat- 
terns are evident, since the values of P correspond to the 
higher-order peaks in Fig. [H (b). The population of the 
m = 1 state, / |^i| 2 dr/7V, is 0.01 in Fig. (a) and 0.05 
in Figure [5] (b). The population of the m = — 1 state is 
very small ~ 10 -4 . 



IV. GAUSSIAN VARIATIONAL ANALYSIS 

To qualitatively examine the Bogoliubov spectra ob- 
tained in Sec. we perform Gaussian variational anal- 
ysis. The variational wave function for the m = state 
has the form 



^ ± ^ eXP V 2d\ 2d* V 



(9) 



where d± and d z are the variational parameters char- 
acterizing the size of the condensate in the radial and 
axial directions. Substituting Eq. (O into Eq. (0) and 



the mean-field energy 

E= J 'drr {-^ 2 + V +!h/> | 2 )vo, (10) 
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where A = cu z /lu±, d± = d±/a±, d z = d z /a±, and 
9o = g N/[(27r) 3 / 2 hu ± al] with a ± = [%/ {Muj^)] 1 ' 2 . 
The variational parameters d± and d z are determined 
so as to minimize Eq. ([T2]) . 

For simplicity, we restrict the magnon excitation to the 
form, 

1>±i(r,t) = e-^l % [a ±1 e-^ t X ±(r) + ^e^x^r)} , 

(13) 

with 
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which corresponds to the lowest mode of L = 1 in Eq. (jSJ). 
Substitution of Eqs. ©, ([TTj). and into Eq. © yields 
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with (ja = gdN/[6(27r) 1 ^ 2 hu;±a 3 L \. Diagonalizing 
Eq. ((TSJ), we obtain the excitation frequency as 

Co 2 = P 2 +A 2 + 2(G + L>i)A 

±2y/ [A 2 + 2{G + L>i)A]P 2 + A 2 D 2 . (20) 

For the parameters of 87 Rb in Fig. [TJ uj in Eq. ([20]) 
becomes imaginary between P ~ 49.9 and 50.01 and the 
maximum value of Im c<; is ~ 0.01. For the parameters of 
23 Na in Fig. [H uj becomes imaginary between P 0.125 
and 0.13 and the maximum value of Im uj is ~ 0.002. 
These results are in qualitative agreement with the first 
peaks of L = 1 in Figs. [1] and [H The differences between 
the variational and numerical results come from the forms 
of the variational wave functions assumed in Eqs. ([9]) and 
(p~3|) : more appropriate variational functions are needed 
for quantitative explanation of the numerical results. 

V. CONCLUSIONS 

In conclusion, we have studied the magnetization dy- 
namics caused by the MDI in a spin-1 BEC in the m = 
hyperfine state prepared in a pancake-shaped trap and 



a magnetic field applied in the axial direction. We 
found that transverse magnetization develops due to the 
MDI breaking the chiral or axial symmetry, and a vari- 
ety of magnetization patterns appear depending on the 
strength of the applied magnetic field. We showed that 
these phenomena occur in spin-1 87 Rb and 23 Na BECs. 
We also performed Bogoliubov analysis and found that 
the initial fluctuations in the magnetization are exponen- 
tially amplified by the dynamical instability. A Gaussian 
variational analysis provided a qualitative explanation of 
the results. 

Our study has shown that magnetization due to the 
MDI strongly depends on the shape of the system. Mag- 
netization dynamics for various trapping potentials in- 
cluding cigar-shaped and lattice potentials merit further 
study. 
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